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We have continued our study of the phase diagram of high temperature QCD with three flavors of improved 
staggered quarks. We are performing simulations with three degenerate quarks with masses less than or equal to 
the strange quark mass and with degenerate up and down quarks with masses mu,d less than the strange quark 
mass. For the quark masses studied to date, we find a crossover that strengthens as m„,d decreases, rather than 
a bona fide phase transition. We present new results for the crossover temperature extrapolated to the physical 
value of mu,d, and for quark number susceptibilities. 



1. THE PHASE DIAGRAM 

The MILC Collaboration is studying high tem- 
perature QCD with three flavors of improved 
staggered quarks 1 using the Asqtad quark ac- 
tion [2]. Simulations are being carried out with 
lattice spacings 1/4T, l/6r and 1/8T. We are 
considering two cases: 1) all three quarks have 
the same mass niq; and 2) the two lightest quarks 
have equal mass niu^d, while the mass of the 
third quark is fixed at approximately that of the 
strange quark m^. We refer to these cases as 
Nf = 3 and Nf ==2 + 1, respectively. We have 
carried out studies with 0.2 nis < fUq < rus for 
the Nf = 3 case, and with 0.1 < rriu^d £ 'm-s 
for Nf = 2 + 1. At the masses we have studied 
to date, we find rapid crossovers, which sharpen 
as the quark mass is reduced, rather than a bona 
fide phase transition. During the past year we 
have extended and clarified our calculations of 
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the quark number susceptibilities, which provide 
an excellent signal for the crossover, and which 
are directly related to event by event fluctuations 
in heavy ion collisions. Comparison of suscep- 
tibilities at different lattice spacings show that 
our results are close to the continuum values. We 
have also sharpened our estimates of the crossover 
temperature for both the Nf — 3 and Nf = 2 + 1 
cases. 

For both the Nf = 3 and 2 + 1 cases we have 
carried out thermodynamics studies on lattices 
with four, six and eight times slices, and aspect 
ratio Ns/Nt — 2. Here Ng and Nt are the spatial 
and temporal dimensions of the lattice in units of 
the lattice spacing. At the lightest quark masses, 
we also performed some simulations with aspect 
ratio three, and obtained results that are indistin- 
guishable from those with aspect ratio two. For 
both cases the standard thermodynamic quanti- 
ties show a crossover from confined behavior at 
low temperature to deconfined behavior at high 
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Figure 1. The real part of the Polyakov loop on 
16^ X 8 lattices for iV/ = 3. 



temperature, as is illustrated in Fig. ^ where we 
plot the real part of the Polyakov loop for Nf — 3 
on 16'^ X 8 lattices. There is a slight trend for 
the temperature dependence of the Polyakov loop 
to be steeper for larger quark masses. This is 
to be expected, since at sufficiently large quark 
masses, it is a bona fide order parameter. The 
ipij^ susceptibility Xtot provides a clear signal for 
the crossover. It is given by 



(1) 



We plot this quantity as a function of temperature 
on 8^ X 4 lattices in Fig. El Note the increase 
in the height of the peak as the quark mass is 
decreased. Our work and that of the Bielefeld 
group |H] strongly suggests that in the real world, 
Nf = 2 + 1, there is no phase transition at the 
physical quark masses. With this assumption, we 
have estimated the critical temperature for Nf — 
2 + 1 at the physical value of rriu^d through an 
extrapolation of the form 



ri Tc = Co + ci [m^/mp)'^ + C2 {aTt,f 



(2) 



where we evaluated Tc for each value of Nt and 
"ni-n/'nip for which we have made measurements 
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Figure 2. The ?/''0 susceptibility as a function of 
temperature on 8'^ x 4 lattices for Nf = 2 + 1. 

from the peak in the susceptibility. For a 
second order phase transition in the 0(4) univer- 
sality class at niu^ = 0, d — 2/pS ^ 1.08. We 
find that Tc = 169(12)(4) MeV with a of 2.1 
for 11 degrees of freedom. The first error is the 
fit error, the second from the uncertainty in ri, 
taken as 0.317(7) fm jlj. To test the sensitivity of 
Tc to d, we have also performed a fit with d = 2, 
which yields = 174(11) (4) MeV with a of 
1.5 for 11 degrees of freedom. So, the goodness of 
the fit does not allow us to prefer either of them. 

2. QUARK NUMBER SUSCEPTIBILI- 
TIES 

In order to study the quark number suscep- 
tibilities j5l6j . we introduce chemical potentials 
fia coupled to a set of mutually commuting con- 
served charges Qa- The quark number suscep- 
tibilities are related to event-by-event fluctua- 
tions in heavy ion collisions |Zj by the fluctuation- 
dissipation theorem 



XcAT) = {{Qa ~ {Qa))iQ0 - (Q/3))) 



(3) 



We work at = 0, so the brackets, (), in Eq. Q 
indicate averages weighted by the standard, real 
Euclidean action for QCD, and (Qa) = 0. 
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Rather than choosing the three independent 
charges to be the number operators for up, down 
and strange quarks, it appears more physical to 
take them to be the z-component of isospin, Qi , 
the hypercharge, Qy, and the baryon number, 
Qb- The rows and columns of the susceptibil- 
ity matrix, Xa,/3(T') are then labeled by /, Y and 
B. In the Nf = 3 case, where niu = rud = nis, 
X is a diagonal matrix, and there are no corre- 
lations among fluctuations in Qj, Qy and Qb, 
while for the Nf = 2-1-1 case, where m„ = rud, 
the only correlations are between fluctuations in 
hypercharge and baryon number. For tempera- 
tures below the phase transition or crossover, the 
lightest particle that can be excited by a chemi- 
cal potential coupled to Qi is the pion, while for 
chemical potentials coupled to Qy and Qb it is 
the kaon and the nucleon, respectively. Above 
the transition temperature each of the chemical 
potentials can excite quark states that are much 
lighter than hadrons, so we expect the diagonal 
elements of x to increase sharply in the vicinity 
of the transition, and they do. This is illustrated 
in Fig. 01 where we plot the diagonal elements of 
the susceptibility matrix, xij, Xy,y and xb,b, 
as a function of temperature for two light quarks 
with mass 0.2 nis and one heavy quark with mass 
nis on 12^ X 6 lattices, xy.y and xb.b have been 
multiplied by factors of 3/4 and 3/2 respectively, 
so that the quantities plotted approach the same 
high temperature limit as xi,i- Also shown is 
Xy,b, the only non-zero off-diagonal matrix ele- 
ment of X for m„ — TTid- It measures correla- 
tions between fluctuations in the hypercharge and 
baryon number. The cocfflcient of xy.b in this fig- 
ure is the geometric mean of those for xy,y and 
Xb.b- The good agreement between these results, 
and similar ones on lattices with eight time slices, 
indicate that they are close to their physical val- 
ues. 
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Figure 3. The non- vanishing elements of the sus- 
ceptibility matrix for two light quarks of mass 
0.2 ms and one heavy quark of mass rris on 12'^ x 6 
lattices. 
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